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BACKGROUND 

Field of the Invention 

[0001] The present invention relates generally to frequency translation. More specifically, the 

present invention relates to frequency translation in a communication system that can be easily 
implemented in digital hardware. The present invention further relates to use of frequency 
f=i translators in numerous applications, including upconverting baseband signals to RF frequencies 
yS for transmission through modems, correcting for local oscillator offsets and FSK modulation. 
If | Background of the Invention 

[0002] Frequency translation is required for numerous applications. For example, transmission 

!L of a baseband signal through a modem requires that the baseband signal be upconverted to a 

J2 suitable transmission frequency. For another example, in wireless modem applications, suitable 

p transmission frequencies are in the RF range. RF synthesizers exist to generate such frequencies. 

However, generally they are tuned to operate at frequencies separated by discrete steps. Due to 
the difficulty of designing RF synthesizers with small step values without severely degrading 
performance due to phase noise (z.e., distortion), RF transmitters use large steps between RF 
frequencies. Thus, it may be difficult to obtain an arbitrary RF frequency using an RF 
synthesizer. 

[0003] Another problem with frequency synthesis arises from the use of local oscillators. Local 

oscillators are often inaccurate. While the amount of inaccuracy can be determined through 
calibration, tuning the oscillator to compensate for the inaccuracy is challenging. Further, if the 
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oscillator is not properly tuned, the application in which the oscillator is being used may not 
function properly. 

[0004] As digital processors become more powerful, functions traditionally performed by analog 

radio frequency (RF) circuitry are migrating to digital hardware and software. One such function 
is frequency translation. Figure 1 is a schematic diagram of a conventional frequency translator. 
A phase accumulator 102 outputs a phase angle 9. Phase accumulator 1 02 is stepped by a phase 
increment A0 in an iterative manner. Thus, phase accumulator 102 increments by phase 
increment A0 on each iteration. The size of the phase increment AO determines the 
corresponding frequency. Phase angle 6 is used as an index into a sin/cos lookup table 104. The 
5 output of sin/cos lookup table is a complex value e 56 . The output of sin/cos table 1 04 provides 
I! one input to a complex multiplier. The other input of the complex multiplier is the signal, X, that 
D is to be translated in frequency. The output of the complex multiplier at each iteration k is 
U y k = e JkAe . X k . Thus, output Y corresponds to input X shifted in frequency by an amount 

S proportional to phase increment A0. Phase accumulator 102 and sin/cos table 104 comprise what 

3 is commonly referred to as a numerically controlled oscillator (NCO) 107. 
f SUMMARY OF THE INVENTION 

[0005] The present invention replaces the numerically controlled oscillator of conventional 

frequency translators with a CORDIC phase rotator. A phase accumulator, coupled to the 
CORDIC phase rotator is initialized with the desired phase angle of rotation. The CORDIC 
phase rotator rotates an input vector by the phase angle of rotation (or a close approximation 
thereto) using only shift and add functions. This enables the present invention to be readily 
implemented in hardware using, for example, an application specific integrated circuit (ASIC). 
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Further, the CORDIC phase rotator can be implemented as a pipeline of stages, each stage 

adding or subtracting a portion of the phase angle of rotation. 
[0006] In addition, the present invention applies stochastic (probabilistic) rounding to the output 

of the CORDIC processor to minimize truncation noise. Stochastic rounding is accomplished by 

adding a 4 bit pseudo random number, generated by a 4-bit pseudo random number generator to 

the least significant 4 bits of the output of CORDIC phase rotator. The carry of that addition is 

added to the 10 most significant bits of the CORDIC phase rotator. 
[0007] Further, when the signal is to be applied to a component that exhibits a nonlinearity such 

as a D/A converter, the present invention applies a dithering technique to the rounded signal to 
1 compensate for the nonlinearity. A dither signal is generated by a 1 -bit random number 

St generator. The signal so generated is filtered by a transversal low pass filter (LPF) to 

q concentrated the energy of the dither signal near zero frequency. The filtered dither signal is 
M added to the signal to be converted to analog. 

|§j)08] The dither signal is removed from the resulting analog signal in two ways. First, the 

O aliasing characteristics inherent in D/A converter operation removes some of the dither signal. 

F 1 Second, band pass filtering removes a substantial portion of any dither signal that remains. 

[0009] The present invention is described in greater detail in the detailed description of the 

invention, the appended drawings and the attached claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0010] Figure 1 is a schematic diagram of a conventional frequency translator. 

[001 1] Figure 2 is a schematic diagram of a frequency translator according to a first embodiment 

of the present invention. 
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[0012] Figure 3 is a schematic diagram of a frequency translator according to a second 

embodiment of the present invention. 
[0013] Figure 4 is a schematic diagram of an exemplary stage of a CORDIC phase rotator. 

[0014] Figure 5A is a frequency domain plot illustrating truncation noise. 

[0015] Figure 5B is a frequency domain plot illustrating truncation noise reduction using 

stochastic rounding. 

[0016] Figure 6 is a schematic diagram of an apparatus for performing stochastic rounding. 

[0017] Figure 7 A is a frequency domain plot illustrating harmonics created due to D/A converter 

nonlinearity. 

jj}18] Figure 7B is a frequency domain plot illustrating reduction of harmonics created by D/A 

fy converter nonlinearities using a dither signal. 

P019] Figure 8 is a schematic diagram of an apparatus for adding a dither signal to a signal to 

compensate for D/A converter nonlinearity. 
|jf020] Figure 9 is a frequency domain plot illustrating a nulling effect inherent in D/A converter 

~2 operation. 

P DETAILED DESCRIPTION OF THE INVENTION 

[0021] The present invention optimizes hardware frequency translation by performing frequency 

translation using only shift and add operations. In addition, the present invention overcomes the 
problems associated with RF synthesizers that have wide steps. The RF synthesizer can be used 
to get close to the desired frequency, and the present invention can add any fine tuning required 
to generate a desired frequency. For example, the present invention can tune in 1 Hz steps, 
which is difficult to accomplish using an RF synthesizer. Similarly, the present invention can be 
used to compensate for offsets in local oscillator frequencies. This is accomplished by 
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calibrating the local oscillator to determine the offset and then using the present invention to 
translate the oscillator-generated frequency to compensate for the offset. 

[0022] The present invention can also be used as a frequency modulator. For example, to 

perform FSK modulation, a carrier frequency applied to the present invention is shifted in 
frequency depending on the value of the applied bit stream. In a similar manner, the present 
invention can be used to perform frequency hopping. 

[0023] Figure 2 is a schematic diagram of a frequency translator 201 according to a first 

embodiment of the present invention. In Figure 2 phase accumulator 102 outputs a phase angle 
to a CORDIC phase rotator 202. CORDIC phase rotator rotates an input X by a phase angle 0. 

O The output Y can be expressed as Y k = e JkA& . X k , which is equivalent to the output of 
| if conventional frequency translators. However, in frequency translator 201 of the present 

Tt invention, sin/cos table 1 04 and complex multiplier 106 are replaced by CORDIC phase rotator 

T 202. 

|[024] The Coordinate Rotation Digital Computer (CORDIC) is a well-known algorithm, which 

-XI was developed in the 1950s, and is described in Andraka, Ray, "A survey of CORDIC 

H algorithms for FPGA based computers," Proceedings of the 1998 ACM/SIGDA 6 th International 

symposium on FPGAs, February 22-24, 1998, Monterey, CA, pp. 191-200. The advantage of 
CORDIC phase rotator 202 is that it uses discrete logic rather than multipliers and lookup tables 
as required for operation of NCOs. Thus, computationally efficient shift and add operations 
replace more computationally intensive multiplies. As a consequence, CORDIC phase rotator 
202 can be implemented in pipelined hardware for maximum speed. 
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[0025] 



[0026] 



p27] 



[0028] 



In essence, the CORDIC phase rotator performs a rotation using a series of specific 
incremental rotation angles selected so that each is performed by a shift and add operation. The 
operation of CORDIC phase rotator 202 is described in more detail below. 

Rotating a complex vector by a phase angle 0 can be accomplished by multiplying the 



complex vector by a rotation matrix A{6) = 



cos# -sin# 
sin# cos# 



Through basic trigonometric 



substitution, A(9) can be represented as A(0) = 



1 



1 -tan0 



. By restricting the 



Vl + tan 2 #L tan # 1 
value of phase angle 9 to arctangent of 2 to some power (i.e., 0 k - tan _1 (2~* ), where k is any 
integer), A(9) can be advantageously written as: 



Vl + 2' 



■2k 



1 

2~ k 



■2~ k 
1 



A(-G k ) = 



■2k 



-2" 



(1) 



Positive or negative rotations can be performed as described below. The 



Vl + 2- 2 * 

value is a scale factor for the CORDIC algorithm. The scale factor can be disregarded for most 
applications. Its only effect is to provide some gain to the CORDIC algorithm. As can be seen, 
using the CORDIC algorithm allows vector rotation using only shifts and adds, provided rotation 
angles are restricted as provide above. Thus, for each k, the input vector X is rotated by an angle 

equal to tan _1 (2 _i ). 

By cascading the rotation matrices, rotations of arbitrary angles (or close approximations 
thereto) can be obtained. The desired angle of rotation can then be obtained by appropriate 
combination of sums and differences of the angles by which each stage can rotate its input. For 
example, to rotate an input vector by 45 degrees, a rotation matrix that rotates the input vector by 
15 degrees can be cascaded with a rotation matrix that rotates the input vector by 30 degrees. 
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Thus, rotation by an arbitrary angle (or a close approximation thereto) is accomplished by the 
relationship A(&)X « A(±0 } )A(±8 2 ) . . . A(±0 n )X . In this relationship successive stages of the 
cascade rotate the vector input to that stage in a positive or negative direction depending upon 
the value of 9 for that stage. 

[0029] Figure 3 is a frequency translator according to a second embodiment of the present 

invention using a 14-stage CORDIC phase rotator. An input vector X is to be rotated by a phase 
angle 0. The inphase (I) and quadrature (Q) components of input vector X are input to a 
CORDIC phase rotator 302. In the preferred embodiment, CORDIC phase rotator 302 comprises 
14 stages, stage 1 through stage 14. Stages 2 through 14 input the inphase and quadrature 

m components of a vector rotated by the previous stage and output a vector rotated in according to 

; M the relation 6 k = tan" 1 (l~ k ), where k is the stage number and 0^ is the amount of rotation 

p. provided by that stage. If 0^ is negative, the angle is rotated in the negative direction for that 

s stage. In addition, each of these stages outputs an adjusted phase angle. The phase angle is 

ID adjusted to account for the rotation provided by that stage. This can be either a positive or 

J 1 negative adjustment depending on the sign of the phase angle input to the stage. For example, if 

p " the desired rotation is 38 degrees and stage 1 rotates the input vector by 25 degrees, then stage 1 

outputs a phase angle of 13 degrees to stage 2. If stage 2, rotates the angle by 1 5 degrees, then 
stage 2 outputs a phase angle of -2 to stage 3. In this example, stage 3 then rotates its input 
vector in the negative direction and adds the amount of rotation to the phase angle. 
[0030] More or fewer stages can be used in CORDIC phase rotator 302 depending on the 

accuracy desired for arbitrary rotation phase angles. Stage 1 inputs the input vector X and the 
user entered or predetermined phase angle 0. Stage 14 inputs the preceding stage's vector and 
rotates it by an amount equal to +/- 0 14 . Stage 14 outputs the input vector X rotated by 
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approximately the phase angle 9. Preferably, each stage uses 14-bit arithmetic for operations on 
the I and Q inputs and 12-bit arithmetic for operations involving phase angle 0, 
[0031] The angle of rotation (phase angle 9) is provided by phase accumulator 304. The angle of 

rotation can be a predetermined angle or user entered. At each stage the input vector is rotated in 
a positive or negative direction so that the total rotation at the output of the 14 th stage is equal to 
or approximately equal to the desired angle of rotation. In addition, the phase angle output by 
each stage is the phase angle input to that stage adjusted to account for the rotation provided by 
that stage. As described above, the process of rotation using the CORDIC algorithm is reduced 
to a series of shifts and adds. 
[J$f>32] Figure 4 is a schematic diagram of an exemplary stage of CORDIC phase rotator 302 to 

nj more clearly illustrate the required shifts and adds. Each stage comprises three adders 402, 404 

D ; and 406. Adder 402 produces the quantity I±2~ k Q. This is the I term provided to the 

g subsequent stage. Adder 404 produces the quantity Q ± 2~ h I . This is the Q term provided to the 

fit subsequent stage. The quantity I +/- jQ output by stage 14 represents the input I +/- jQ rotated 

yp by phase angle 9 (or a close approximation thereto). 

f0O33] It would be apparent to those skilled in the art that the I and Q terms multiplied by the 2" k 

is equivalent to a right shift of the respective I or Q term by k places. As described above, k is 
equal to the stage number. The shift is represented in Figure 4 by shift operators 408 and 410. 
Shifting operations can be performed efficiently in digital hardware. 

[0034] Adder 406 outputs the phase angle to input to the next stage. This phase angle is the sum 

of the input phase angle and the angle supplied by that stage (?>., the phase angle adjustment 
described above). The output of adder 406 is minimized. Thus, whether to add or subtract a 
stage's phase angle is controlled by the most significant bit (MSB) of the phase angle input to 



-8- 



that stage. As is well known to those skilled in the art, the MSB indicates whether the value is 
positive (corresponding to an MSB of 0) or negative (corresponding to an MSB of 1). 
Consequently, if the phase angle input to a stage is positive, the rotation angle provided by that 
stage is subtracted from the phase angle input that stage. Conversely, if the phase angle input to 
a stage is negative, the rotation angle provided by that stage is added to the phase angle input that 
stage. The output of adder 406, 0±<f> k , (where 0 is the phase angle input to a stage and <p h is the 
angle that the k* stage corresponds to) is passed to the next stage as the next stage's input phase 
angle. 

[0035] The MSB of the phase angle input to a particular stage is used to control whether adders 

yp 402 and 404 add or subtract their inputs. This follows from the relationships represented by 

fll equation (1) above. Inverter 412 is used to maintain the relationship of equation (1). 

0H)36] One concern with using CORDIC is spurious outputs. If 10 bits are used to represent the 

f in-phase and quadrature portions of the input signal, then it is undesirable to use 10-bit arithmetic 

m operations to perform the shift and adds described above. One reason is that simply shifting 

=|j values to the right is effectively a truncation of the rightmost bits of the value. Truncation adds 

error that can be avoided. As is well known, this error manifests itself as spurious noise in the 
spectrum of the output of a process employing truncation. For example, Figure 5A illustrates an 
exemplary spectrum of a 500 Hz sine wave generated by truncating to 10 bits. Spike 502 
corresponds to the power in the sine wave at 500 Hz. However, there is also a significant spike 
504 observed at 2500 Hz. Spike 504 represents truncation noise. 
[0037] The present invention removes the contribution of spurious noise by using more precision 

to perform the arithmetic required for the CORDIC phase rotator than the input precision. For 
example, in the present invention, the input signal has a precision of 10 bits, but 14 bits of 
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precision are used in the arithmetic required for the CORDIC phase rotator. This technique 

reduces truncation error for 10 bits of precision if ensuing operations retain the higher precision, 

Le. 9 continue to use 14 bits of precision. 
[0038] However, not all applications retain such precision. For example, if a 10-bit D/A 

converter follows the CORDIC phase rotator, the 14 bits must be truncated to 10 bits. This 

results in truncation error as described above. 
[0039] Figure 6 is a schematic diagram of an apparatus for reducing the error due to truncation 

through randomization of the error using stochastic rounding. The apparatus illustrated in Figure 

6 can also be implemented as a method as would be apparent to those skilled in the art. In the 
yg present case, truncation from 14 bits to 10 bits is required, but the technique described can be 

f ll extended as required for a particular implementation. A 4-bit random number generator (RNG) 

Cl 602 generates a 4-bit random number. The 4-bit random number is added to the 4 least 

^ significant bits (LSBs) of the output of the CORDIC phase rotator in adder 604. The carry bit is 

J;! added to the 1 0 MSBs of the 14-bit number output by the CORDIC phase rotator in adder 606. 

,/i The result of the addition performed by adder 606 is taken as the 10-bit output value of the 

p CORDIC phase rotator. 

[0040] The rounding apparatus illustrated in Figure 6 performs a stochastic (probabilistic) 

rounding function. Thus, the expected value of the output of the function is equal to the input. 
That is, the average error at the output of the function is 0. Thus, while truncation noise has not 
been completely eliminated, it has been significantly reduced, and spread out, rather than being 
concentrated. 
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[0041] Figure 5B illustrates reduction of spurious noise by using the apparatus described in 

Figure 6. In this case, a spurious peak 508 at approximately 2500Hz is significantly reduced as 
compared to the corresponding spurious peak 504 in figure 5A, where only truncation is used. 

[0042] The signal that has been generated at the output of the CORDIC phase rotator is 

translated in frequency. It can now be input into a D/A converter, for example, to generate a 
signal suitable for transmission by a modem. However, D/A converters are not ideal. One 
aspect of the non-ideal nature of D/A converters is a non-linear response to an input. 
Consequently, if a sine wave is input to a D/A converter, it will include undesired harmonics of 
the sine wave frequency. Figure 7A is an exemplary spectrum of an output of a D/A converter 

.ft when a sine wave is input to it. The sine wave has a frequency of approximately 1 300 Hz as 

fy shown by the largest peak 702. However, undesirable harmonics 704a, 704b and 704c are 

Q clearly evident. 

(0043] Figure 8 illustrates a schematic diagram of an apparatus for largely eliminating the 

harmonics caused by the D/A nonlinearity typically present in most D/A converters. The 
,fi apparatus of Figure 8 can also be implemented as a method as would be apparent to those skilled 

y* in the art. Using an adder 802, a random signal is added to the rounded signal data generated by 

the stochastic rounding process described above. The random signal is generated by a 1-bit 
random number generator (RNG) 806, The signal is bandlimited by filter 804 to a part of the 
frequency spectrum that will not be used. This allows the added random signal to be filtered out 
later. The technique illustrated in Figure 8 is known as dithering. Like the rounding technique 
described above, dithering does not completely eliminate the error due to the nonlinearity, but 
rather, randomizes it. 
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[0044] LPF 804 is preferably a transversal filter. As is well-know, the output of a transversal 

N 

filter is given by the relation, Y k =^ K^i-k > where h^ is the coefficient of the filter tap, and 

Jfc=0 

Xi-k represents each successive bit of the 1-bit pseudo-random sequence generated by RNG 806, 
Using well-known techniques, the coefficients of the filter taps are chosen to bandlimit the 
broadband output of RNG 806 to a part of the frequency spectrum that is not used. As explained 
above, this allows the random (dither) signal to be filtered out later. 

[0045] Filter 806 can be any filter. However, it is optimal to make filter 806 a low pass filter 

(LPF). This allows the random signal to be placed near zero frequency (/.<?., DC). This is 

CI optimal because it enables the aliasing effect of a D/A converter, such as D/A converter 808, to 

Hi effectively null out the random signal The desired signal is optimally placed at 1/4 of the 

%i sample rate. 

jQ046] Figure 7B illustrates the frequency spectrum that results from the operation of the 

O apparatus shown in Figure 8. The spurious noise present at the harmonics has been effectively 

□ spread out over a large bandwidth. The desired signal energy is peak 702. The dither signal 

H energy is distributed across the hump 706. The small noise humps 708a-n result from the 

sidebands of LPF 804, Choosing appropriate coefficients for LPF 804 can make hump 706 
narrower and/or reduce the power in the sidebands 708a-n. For reasons described below, 
operation of LPF 804 centers the dither energy at approximately zero frequency. 
[0047] Figure 9 illustrates the nulling effect of the D/A converter. The desired signal 902a is 

placed a 1/4 of the sample rate. Figure 9 illustrates the sinx/x response of a typical D/A 
converter. Aliasing, inherent in D/A converter operation, causes images of the desired signal to 
be created. These are images 902b, 902c, 902d, 902e, 902f, 902g and 902h. As can be seen, 
placing desired signal 902 at approximately 1/4 the sample rate places it near a peak in the sinx/x 
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function of the D/A converter in each of the ensuing images. Thus, the image is not significantly 
attenuated, 

[0048] The dither signal 904 is placed close to zero by operation of LPF 806. D/A converter 

aliasing also causes images of the dither signal to be created. These images are peaks 904b, 
904c, 904d and 904e, It is apparent however, that the images of the dither signal appear at or 
very close to the nulls of the D/A converter's frequency response. Consequently, the dither 
signal is largely eliminated by the inherent nature of the D/A converter. As shown by hump 706 
in Figure 7B, the dither signal energy is wideband with a width depending on the cutoff 
frequency of filter 804. Consequently, some of its energy is likely to be present though severely 

if| attenuated by the sinx/x response of the D/A converter. 

P049] In modem applications, one of the images of the D/A converter is typically chosen for 

p transmission. For example, if the sample rate is 18.432 MHz, the desired signal is placed at 

^ substantially 4.608 MHz (i.e., 1/4 the sample rate). By operation of the D/A converter, an image 

5[ suitable for transmission is present at 32.256 MHz by operation of the D/A converter. The 

J corresponding dither signal image is at approximately 36.864 MHz, which is a null of the D/A 

M converter frequency response for an 1 8.432 MHz sample rate. Consequently, the dither signal is 

largely eliminated in the ensuing transmission. A suitable analog bandpass filter can be designed 
to pass the desired signal after nulling by the D/A. For example, the image 902d can be used and 
the corresponding image of the dither signal 904c is largely suppressed by operation of the D/A 
converter. 

[0050] The foregoing disclosure of the preferred embodiments of the present invention has been 

presented for purposes of illustration and description. It is not intended to be exhaustive or to 
limit the invention to the precise forms disclosed. Many variations and modifications of the 

-13- 



embodiments described herein will be obvious to one of ordinary skill in the art in light of the 
above disclosure. The scope of the invention is to be defined only by the claims appended 
hereto, and by their equivalents. 

Further, in describing representative embodiments of the present invention, the 
specification may have presented the method and/or process of the present invention as a 
particular sequence of steps. However, to the extent that the method or process does not rely on 
the particular order of steps set forth herein, the method or process should not be limited to the 
particular sequence of steps described. As one of ordinary skill in the art would appreciate, other 
sequences of steps may be possible. Therefore, the particular order of the steps set forth in the 
specification should not be construed as limitations on the claims. In addition, the claims 
directed to the method and/or process of the present invention should not be limited to the 
performance of their steps in the order written, and one skilled in the art can readily appreciate 
that the sequences may be varied and still remain within the spirit and scope of the present 
invention. 
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